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Abstract. Inspired by the study of community structure in connection net- 
works, we introduce the graph polynomial Q (G; y), the bivariate generating 
function which counts the number of connected components in induced sub- 
graphs. 

We give a recursive definition of Q {G; x,y) using vertex deletion, vertex 
contraction and deletion of a vertex together with its neighborhood and prove 
a universality property. We relate Q (G; x, y) to other known graph invariants 
and graph polynomials, among them partition functions, the Tutte polynomial, 
the independence and matching polynomials, and the universal edge elimina- 
tion polynomial introduced by I. Averbouch, B. Godlin and J.A. Makowsky 
(2008). 

We show that Q (G; x, y) is vertex reconstructiblc in the sense of Kelly and 
Ulam, discuss its use in computing residual connectedness reliability. Finally 
we show that the computation of Q (G; y) is jJP-hard, but Fixed Parameter 
Tractable for graphs of bounded tree- width and clique-width. 
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1. Introduction 

1.1. Motivation: Community Structure in Networks. In the last decade 
stochastic social networks have been analyzed inatheniatically from various points 
of view. Understanding such networks sheds light on many questions arising in 
biology, epidemology, sociology and large computer networks. Researchers have 
concentrated particularly on a few properties that seem to be common to many 
networks: the small-world property, power-law degree distributions, and network 
transitivity. For a broad view on the structure and dynamics of networks, see |36) . 
M. Girvan and M.E.J. Newman, [21], highlight another property that is found 
in many networks, the property of community structure, in which network nodes 
are joined together in tightly knit groups, between which there are only looser 
connections. 

Motivated by [35], and the first author's involvement in a project studying social 
networks, we were led to study the graph parameter (G), the number of vertex 
subsets X C V with i vertices such that G [X] has exactly j components, qij (G), 
counts the number of degenerated communities which consist of i members, and 
which split into j isolated subcommunities. 

The ordinary bivariate generating function associated with (G) is the two- 
variable graph polynomial 

n n 

Q(G;x,y) = ^5]g,y (G)xV^. 

1=0 j=Q 

We call Q (G; x, y) the subgraph component polynomial of G. The coefficient of y'^ 
in Q (G; x, y) is the ordinary generating function for the number of vertex sets that 
induce a subgraph of G with exactly k components. 

1.2. Q{G;x,y) as a Graph Polynomial. There is an abundance of graph poly- 
nomials studied in the literature, and slowly there is a framework emerging, [301 
1311 125] , which allows to compare graph polynomials with respect to their ability to 
distinguish graphs, to encode other graph polynomials or numeric graph invariants, 
and their computational complexity. In this paper we study the subgraph compo- 
nent polynomial Q (G; x, y) as a graph polynomial in its own right and explore its 
properties within this emerging framework. 

Like the bivariate Tutte polynomial, see \\^ Chapter 10], the polynomial Q (G; x, y) 
has several remarkable properties. However, its distinguishing power is quite dif- 
ferent from the Tutte polynomial and other well studied polynomials. 

Our main findings are: 

• Q (G; X, y) distinguishes graphs which cannot be distinguished by the match- 
ing polynomial, the Tutte polynomial, the characteristic polynomial, or the 
bivariate chromatic polynomial introduced in [TB] (Section [3]). 

• Nevertheless, we construct an infinite family of pairs of graphs which cannot 
be pairwise distinguished by Q {G;x,y) (Proposition [20]) . 

• The Tutte polynomial, satisfies a linear recurrence relation with respect 
to edge deletion and edge contraction, and is universal in this respect. 
Q (G; X, y) also satisfies a linear recurrence relation, but with respect to 
three kinds of vertex elimination, and is universal in this respect. (Theo- 
rems [13| and [21]) . 



2 



P. TITTMANN, I. AVERBOUCH, AND J. A. MAKOWSKY 



• A graph polynomial in three indeterminates, ^(G; a;, ?/, z), which satisfies 
a linear recurrence relation with respect to three kinds of edge elimation, 
and which is universal in this respect, was introduced in [HIS]. It subsumes 
both the Tutte polynomial and the matching polynomial. For a line graph 
L{G) of a graph G, we have Q {L{G); x,y) is a substitution instance of 
C(G;x,?/,z) (Theorem [221). 

• For fixed positive integer n the univariate polynomial Q{G;x,n) can be 
interpreted as counting weighted homomorphisms, [50], and is related to 
the Widom-Rowlinson model for n particles (Theorem [TU]) . 

• Q{G;x,y) is reconstructible from its vertex deletion deck in the sense of 
[I21[II] (Theorem [27|). 

• Q{G; X, y) can be used (Section[8|), to compute the probability Pk (G) that a 
vertex induced subgraph of G has exactly k components from the subgraph 
polynomial. For fc = 1 this is known as the residual connectedness reliability 
(Section [8]). 

• Also like for the Tutte polynomial, cf. [27], Q{G;xo,yo) has the Difficult 
Point Property^ i.e. it is (jP-hard to compute for all fixed values of (xq, yo) G 
M.'^ —E where E' is a semi- algebraic set of lower dimension f Theorem [29]). In 
[5T] it is conjectured that the Difficult Point Property holds for a wide class 
of graph polynomials, the graph polynomials definable in Monadic Second 
Order Logic. The conjecture has been verified only for special cases, [6l[7l[8]. 

• Q{G] xq, yo) is fixed parameter tractable in the sense of [19] when restricted 
to graphs classes of bounded tree- width (Proposition [31]) or even to classes 
of bounded clique- width (Proposition [32]). For the Tutte polynomial, this 
is known only for graph classes of bounded tree- width, [371 ^2, j33J . 

Outline of the paper. The paper is organized as follows: In Section [5] we in- 
troduce the polynomial Q{G] x, y) and its univariate versions. In Seection [3] we 
discuss the distinguishing power of Q{G;x,y) and compare this to other graph 
polynomials. In Section [4] we show how certain graph parameters are definable 
using Q(G;x,y) and relate it to partition functions and counting weighted homo- 
morphisms. In Section [S] we give a recursive defintion of Q{G; x, y) using deletion, 
contraction and extraction of vertices and show that Q{G;x,y) is universal. We 
also compare it to the universal edge elimination polynomial ^{G]x,y, z) defined 
in [4l[5] and give a subset expansion formula for Q{G; x,y). In Section [6] we prove 
decomposition formulas for Q{G]x,y) for clique separators. In Section [7] we show 
the reconstructibility of Q{G; x,y). In Section [5] we discuss its use to compute 
the residual connectedness reliability. In Section [9] we discuss the complexity of 
computing Q(G; x, y). In Section [10] we draw conclusions and state open problems. 

2. The Subgraph Component Polynomial Q{G; x, y) 

2.1. The Bivariate Polynomial. Let G — {V,E) be a finite undirected graph 
with n vertices and let fc < n be a positive integer. Assume the vertices of G 
fail stochastic independently with a given probability q = 1 — p. What is the 
probability that a subgraph of G with exactly k components survives? The solution 
of this problem leads to the enumeration of vertex induced subgraphs of G with 
k components. For a given vertex subset X Q V, let G [X] be the vertex induced 
subgraph of G with vertex set X and all edges of G that have both end vertices in 
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X. We denote by k (G) the number of components of G. Let qij (G) be the number 
of vertex subsets X C V with i vertices such that G [X] has exactly j components: 

(7y- (G) ^\{X CV:\X\^tAk{G [X]) = j}\ 

The ordinary generating function for these numbers is the two- variable polynomial 

71 n 

Q(G;a:,2/)=^^g.,(G)xV^'. 

We cah Q (G; x, y) the subgraph component polynomial of G. Since loops or parallel 
edges do not contribute to connectedness properties of a graph, we assume in this 
paper that all graphs are simple. 




Figure 1. The star Stars = ^1,3 

The star ii'1.3, presented in Figure[Tl has the subgraph polynomial 

Q (Ki^s; x,y) ^1 + ixy + Sx'^y + ix^y + x'^y + Sx^y^ + x^y^. 

The term Zx^y"^ tell us that there are 3 possibilities to select two vertices of G that 
are non-adjacent. 

The empty set induces the null graph N — (0, 0) that we consider as being 
connected, which gives goo (G) = Q (G; 0, 0) = 1 for any graph. Substitution of 1 
for y results in an univariate polynomial that is the ordinary generating function 
for all subsets of V, i.e. Q(G; x, 1) = (1 + x)". 

2.2. Univariate Polynomials. The coefficient of y^ in Q (G; x, y) is the ordinary 
generating function for the number of vertex sets that induce a subgraph of G with 
exactly k components: 

(G;x)= [2/'=]Q(G;x,y) 

We call the polynomial Qk for fc G N again subgraph component polynomial. The 
subscript as well as the number of variables should avoid confusion with the formerly 
defined subgraph polynomial. The subgraph polynomial Qi (G; x) is of special 
interest. We rename this polynomial to 

n 

5(G;x) =gi(G;x) =^s,(G) x\ 

It counts the connected vertex induced subgraphs of G. A separating vertex set of 
a connected graph G — {V, E) is a subset X (ZV such that G — X is a disconnected 
graph. 

Theorem 1. Let G = {V,E) be a connected graph with n vertices. Let Ck (G) be 
the number of separating vertex sets of cardinality k for k — 0,1, ...,n. Then the 
coefficients of the subgraph polynomial S (G; x) are given by 



Sk (G) = 
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Proof. If X is a separating vertex set then V \ X induces a disconnected graph. 
Conversely, if X C y is not a separating vertex set of G then G \ X] is connected. 

□ 

We conclude that 2" — S{G; 1) is the number of all separating vertex sets of G. 




Figure 2. Non-isomorphic trees with the same subgraph polynomial 

A graph invariant / is trivial on a class of graphs K of for any two graphs Gi 
and G2 with the same number of vertices we have /(Gi) = /(G2). 

Proposition 2. All non-isomorphic trees with up to nine vertices have different 
subgraph component polynomials. In other words we have: The graph polynomials 
Qk{G;x) and Q{G;x,y) are not trivial on trees. 

However, we have: 

Proposition 3. There exist a unique pair of non-isomorphic trees with 10 vertices 
sharing the same subgraph component polynomial. 

Proof. Figure [5] shows these trees. This statement is true for S (G; x) as well as for 
the (general) subgraph polynomial Q (G; x,y). □ 

In Section Owe shall see how to use this to generate infinite families of pairs of 
graphs which are not distinguished by Q(G; x,y). 

3. Distinctive Power 

We denote by m{G;x) = mi(G)a;* be the matching polynomial with mi{G) 
the number of i-matchings of G, by p(G; x) be the characteristic polynomial, by 
T(G; X, y) the Tutte polynomial, and by F(G; x, y) the bivariate chromatic polyno- 
mial introduced in |18j . 




Figure 3. The graphs Gi, G2 
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Figure 4. The graphs G3, G4 





Figure 5. The graphs G5, Ge 



Proposition 4. For the graphs Gi\i — 1, . . .6 from Figures\^ [7| [3 and 
for P4 and Ki^^ we have 

(1) p(Gi;x) -p(G2;x) but Q{Gi;x.,y) ^ Q{G2;x,y). 

(2) TO(G3;a;) = m(G4;x) but Q{Gz;x,y) ^ Q{Gi;x,y). 

(3) P(G5;x,y) = P(G6;a;,y) g(G5; ^ Q(G6; x, y). 

(4) r(F4;x,y) = r(ifi,3;a:,y) butQ{Pr,x,y)^Q{Kix,x,y). 

Proof. (1) and (2) are easy to verify. 

For (3) P(G5;x,2/) = P(G6;a;,y) is from 1^. For Q(G5;x,y) ^ Q(G6;x,y) we 
compare [x'^y^]Q{G5, x,y) with [a;42/3](5(G6, y). 

For (4) we use that the Tuttc polynomial does not distinguish trees of the same 
size, but that Q{G]x,y) distinguishes all trees of size up to nine vertices, [3 D 



Remark 5. Q{G;x,y) does not distinguish between graphs which differ only by 
the multiplicity of their edges, whereas for the Tuttc polynomial this is not the case. 
Let K^r^^ denote the complete graph with m edges between any two distinct vertices. 
Then we have T{K^; x, y) ^ T{Kl; x, y) but Q{K^; x, y) = Q{K^; x, y). 

Problem 6. Are there simple graphs distinguished by p{G;x), m{G;x), P{G]x,y) 
or T{G; x, y) which are not distinguished by Q{G; x, y) ? 

We say that a simple graph G is determined by a graph polynomial f if for 
every simple graph G' such that /(G) — f{G') we have that G is isomorphic to 
G'. The class of simple graphs K is determined by a graph polynomial f if every 
graph G G if is determined by /. This notion has been studied in [351 HHj for 
the chromatic polynomial, the Tutte polynomial and the matching polynomial. It 
is shown, e.g., that several well-known families of graphs are determined by their 
Tutte polynomial, among them the class of wheels, squares of cycles, complete 
multipartite graphs, ladders, Mobius ladders, and hypercubes. 
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It follows from Proposition [8] that the class of empty graphs En is determined 
by Q{G; x,y, z), and so is the class of complete graphs Kn- Note that, since 
T{En;x,y) — 1 for all n G N, the class of empty graphs is not determined by 
the the Tutte polynomial. It follows from Proposition [3] that the class of trees is 
not determined by Q{G; x, y, z). 

Proposition 7. The class of graphs of the form Starn = K\^n is determined by 
Q{G;x,y,z). 

Proof. It is easy to verify, that if [x^'-^^y]Q{G;x,y) — 1, [x'^y'^]Q{G; x,y) = 1 and 
[a;"+^](5(G; x, y) = 0, then G is isomorphic to Starn- □ 

4. Combinatorial Interpretations 

4.1. Evaluations and Coefficients of Q{G; x,y). For a polynomial f{x,y), let 
[x^y^] f {x, y) be the coefficient of x^y^ in / (x, y) and let deg^, / be the degree with 
respect to the variable x. 

Proposition 8. The following graph properties can be easily obtained from the 
subgraph polynomial: 

(1) The number of vertices: 

n{G) = \V{G)\= deg, Q (G; x, y) = log^ Q (G; 1, 1) 

(2) The number of edges: 

e{G)^[x^y]Q{G-x,y) 

(3) The number of components: 

fc(G)=deg([x"(G)]Q(G;a;,y)) 

Theorem 9. The degree of the subgraph component polynomial Q (G; x, y) with 
respect to y is the cardinality of a maximum independent set ofG (the independence 
number): 

degy Q (G; x,y) = a (G) 

Proof. Let X C be a maximum independent set of G. In this case, we have 
k (G [X]) = \X\ and hence deg^ Q (G; x,y) > a (G). Assume that there exists a set 
Y C V with fc (G [y]) > \X\. Then we obtain an independent set X' by selecting 
one vertex of each component of G [Y] such that \X'\ > |X| ~ a contradiction. □ 

Let Oi (G) be the number of independent vertex sets of size i of G. The indepen- 
dence polynomial of G is defined by 

n 

I{G■,x)=J2a^ (G) x\ 

i=0 

As a consequence of Theorem ^ we can derive the independence polynomial of G 
from the subgraph component polynomial. Let [?/■'] Q (G; x, y) denote the coeffi- 
cient of y^ in Q (G; x, y). This coefficient is a polynomial in x where the coefficient 
of X* counts the vertex subsets of cardinality i of G that induce a subgraph with 
j components. A vertex set X C V is independent if and only if k (G [X]) = \X\. 
Hence [x^] [y^] Q (G; x, y) = Oj (G) is the number of independent vertex sets of size 
3 of G. 
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4.2. Partition Functions. In this subsection we show that the subgraph poly- 
nomial Q{G; x,yo) for any x G M and fixed yo G N can be viewed as a partition 
function, using counting of weighted graph honiomorphisms. Partition functions 
were first studied in the context of statistical physics and have recently attracted 
much attention, [34j . A systematic study of the question which graph invariants 
can be presented as partition functions has been initiated in l23j . 

A weighted graph {H, a, (3) consists of a graph H ~ {V{H), E{H)) with a assign- 
ing weights to vertices and (3 to edges. The partition function ZH.a,i3{G) associated 
with {H, a, (3) is defined by 

h : V ^ Vh v£V {u,v)eE 
homomorphism 




Figure 6. Auxiliary graph Stavy for y — 6 



Let {Stavy, a,/?) be a weighted star with y + 1 vertices and with all loops. The 
central vertex is vq. An example of Stavy for j/ = 6 is shown on Figure [H The 
weight functions are defined as follows: 



1 if V = Vq 

X otherwise 



(i{u,v) = l 

Theorem 10. Let Z(^fj.a,i3){G) he the partition function associated with H = Stary 
and a,P as above. Then, for all nonnegative integers y and all x £ R, we have 

Q{G]X,y) = Z(^stary,a,l3) 

Proof. Let us start with the definition of Zh{G). Under every mapping h : V <—> Vh, 
let A C F be the subset of vertices that are not mapped to vq- Let us count the 
homomorphisms that map the subset A into vi,. . . ,Vy: there are exactly y'^^'^I^D 
such homomorphisms, because every connected component of G[A] must be mapped 
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into a single vertex. Finally, we get 



h -.V v^V 
homomorphism 



= E = E E 

/i : V ^ Vff ACV h : V ^ Vh 

homomorphism homomorphism 

V ^ A '—f h{v) vq 

ACV 

which by ([5]) completes the proof □ 

It is open whether there are other points in which Q(G;x,y) is definable as a 
partition function. 

Remark 11. The auxiliary graph H ~ Starn is called in physical literature The 
Widom-Rowlinson model for n particles. The homomorphisms to Starn are called 
Widom-Rowlinson configurations, [5D] . 

5. Recursive Definition and Subset Expansion 

5.1. Recurrence Relation for Vertex Elimination. We turn now our atten- 
tion to the investigation of properties of the subgraph polynomial that support 
its computation. The first statement concerns the multiplicativity with respect to 
components of the graph. 

Theorem 12 (Multiplicativity). (1) Let G = GiUGi he the disjoint union of 
the graphs Gi and G2. Then 

Q{G; X, y) = Q{Gi;x, y) ■ Q{G2;x, y). 

(2) In particular, if G = {V, E) consists of c components Gi,...,Gc then the 
subgraph polynomial satisfies 

c 

Q{G-x,y)^\{Q{G,-x,y). 

i=i 

Proof. In case c — 2, each subset X QV oi cardinality k is the disjoint union of two 
subsets C F (Gi) and C 1/ (G2) with \Xi\ = j and \X2\ = i~j. The number 
of components of G [X] — G [Xi U X2] is the sum of the number of components of 
G [Xi] and G [X2]. We obtain 

i k 

(1) (G) = E E (Gl) (G2) . 

j=0 1=0 

Thus for a graph with two components, the subgraph polynomial satisfies 

Q{G;x,y) = Q {Gi;x,y)Q {G2;x,y) . 

We obtain the statement of the theorem by induction on the number of components. 

□ 

We distinguish three types of vertex elimination: 



ENUMERATION OF VERTEX INDUCED SUBGRAPHS 



9 




Figure 7. Vertex contraction 

Vertex deletion: For a given vertex v & V (G), let G — v the graph obtained 
from G by removal of v and all edges that are incident to v. We call this 
operation vertex deletion. 

Vertex extraction: Similarly, let G — X be the graph obtained from G by 
removal of all vertices of the set X C V. Let N (v) be the set of vertices 
that are adjacent to w in G (the neighborhood of v). We denote by N [v] the 
closed neighborhood of a vertex v in G, i.e. the set of all vertices adjacent 
to V including v itself. The operation G — N[v] is called vertex extraction. 

Vertex contraction: A further special graph operation is needed here - the 
contraction of a vertex. That is the graph G/v obtained from G by removal 
of V and insertion of edges between all pairs of non-adjacent neighbor ver- 
tices of V. Figure [7] shows an example graph and the graph obtained by 
vertex contraction. 

Theorem 13. Let G = (V, E) be a graph and v Cz V . Then the subgraph polynomial 
satisfies the decomposition formula 

Q (G; x,y) Q {G - v;x,y) + X {y - 1)Q {G - N [v] ■,x,y)+ xQ {G/v; x, y) . 

Proof. Let us first consider all vertex induced subgraphs of G that do not contain 
vertex v. These subgraphs are also vertex induced subgraphs of G — v. Conse- 
quently, 

Q{G~v;x,y) 

enumerates all induced subgraphs not including the vertex v. 

In a second step we count all vertex induced subgraphs that contain vertex 
V but none of its neighbors in G. In this case, the vertex v forms a connected 
component consisting of v only. The rest of the induced subgraph is a subgraph of of 
G — N[v]. All these subgraphs are enumerated by Q (G — [v] ;x,y). However, the 
component built by v contributes one vertex and one component to the polynomial. 
Thus we obtain the generating function 

xyQ {G-N [v] ■,x,y). 

In our enumeration so far we missed exactly those vertex induced subgraphs that 
contain v and at least one of its neighbors together in one component. We include v 
in the corresponding candidate set, remove it from G, and multiply the generating 
function by x (not by xy because we do not increase the number of components). In 
order to trace the components, we have to simulate the paths using v in G. These 
paths are no longer present in G — v. This task is best performed by using G/v 
instead oi G — v. Thus we obtain the contribution xQ {G/v; x, y) to the generating 
function. Unfortunately, this polynomial enumerates induced subgraphs that do 
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not contain any vertices from N (v), too. We can fix this problem by subtraction 
of xQ {G — N [v] ;x,y), which gives 

xQ [G/v; X, y) - xQ{G - N [v] ; x, y) 

as final contribution to the generating function. □ 

Corollary 14. Let v d V be a vertex of degree 1 in G — {V, E) and let w be its 
only neighbor in G. Then 

Q (G; X, y) = {l + x)Q{G~ v; x,y) + x {y ~ 1) Q {G - {v, w} ■,x,y). 

Proof. Notice that in this case G/v — G — v and G — N[v] = G — {v, w}. Then the 
statement follows immediately from Theorem 1131 □ 

5.2. Some Easy Computations. The subgraph polynomial can be easily com- 
puted for certain special graphs. 

Proposition 15. For the complete graphs Kn and the empty graphs En (the com- 
plement of Kn) we have: 

(1) Q{Kn;x,y) = y[l + xf -y+l. 

(2) Q{En-,x,y) - {l + xyf. 

Proof. (1) In a complete graph Kn each vertex subset except the empty set induces 
a connected subgraph. 

(2) In the empty graph £"„ each subset X C V induces a subgraph with \X\ 
components. □ 

From Theorem 1131 we obtain a recurrence relation for the subgraph polynomial 
of the paths P„ : 

Proposition 16. 

Q{Pn;x,y) = {l + x)Q{Pn-i;x,y) + x{y- l)Q {Pn-2\x,y) 

Together with the initial values, 

Q{Po'iX,y) = 1, 
Q{Pi;x,y) = l + xy, 



equation (|16p determines the subgraph polynomial of P„ uniquely. The explicit 
solution is 

Q ^p^.^^y) = + « ( 2x{l-y) y+' 1-x-a px(l-y) ^ 
2a \l + x — aj 2a \l-|-a;-|-a 



with a — ^ I — 2x + x'^ + 4xy. 

The subgraph polynomial of the cycle C„ satisfies another recurrence relation: 

Proposition 17. 

Q {Gn;x,y) = Q {Pn-i;x,y) + x{y-l)Q {Pn^z],x,y) + xQ {Cn-i;x,y) 

The join GV H oi two graphs G = {V, E) and H [W, F) with ^ n 1^ = is the 
graph obtained from G U if by introducing edges form each vertex of G to each 
vertex of H. Consequently, the join of two empty graphs Kg and Kt is the complete 
bipartite graph Ks^t- 
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Theorem 18. Let G = {V,E) and H{W,F) be two graphs with V n W 
\V\ = s, \W\ = t. Then 



'{GVH-x,y)^Q (G; x,y) + Q (H; x, y) + [(1 + xf - 1] (1 + xf - 1 



Proof. All vertex subsets oiV IJW belong to exactly one of three classes: 

(1) subsets of V , 

(2) subsets of W, 

(3) subsets that have at least one vertex of V and at least one vertex of W . 

The first two classes are counted by Q{G;x,y) and Q {H;x,y), respectively. 
The empty set is counted twice, which is corrected by subtracting one. All vertex 
subsets of the third class induce connected subgraphs of G W H. The generating 



function for the number of subsets of this class is [(1 + xY — 1] (1 + x)* — 1 



□ 



From Theorem[THJ we deduce the subgraph polynomial of the complete bipartite 
graph: 



Corollary 19. 

Q (Ks^tix, y) = il + xyY + (1 + xy)' + [(1 + xf - 1] (1 + xf - 1 



y-1 



Propositions [TBI and [T71 and Theorem [T51 are explicit instances of general results 
which follow from the fact that Q{G;x,y) is definable in Monadic Second Order 
Logic. We shall discuss this feature in Section [5T5l 

We can use Theorem [T8l and the multiplicativity of Q{G;x,y) to prove the fol- 
lowing: 

Proposition 20. There are infinite families of pairs of non-isomorphic graphs with 
a fixed number of connected components which are not distinguished by Q{G; x, y). 

Proof. Let G be a graph. We define inductively 

Fo{G) = G 
Jo(G) = G 
F,,+i{G) = F„(G)UG 
J„+i(G) = J„(G)VG 

Let Tri and Tr2 be the two trees from Figure [21 Then, using Theorem [TBI and 
the multiplicativity of Q{G; x, y) we have for all n e N 

Q{Fr,{Tri)-x,y) = Q(F„(Tr2); y) 

and 

Q( J„(Tri); x, y) = Q{Jn{Tr2); x, y) 
For G connected, the graphs Jn{G) are connected. The graphs Fn{G) have exactly 
n components. So for m components we combine Fm-i{G) U Jn{G) which has m 
components. □ 

5.3. The UniversaUty Property of Q {G; x,y). The vertex decomposition for- 
mula represented in Theorem [T3| can be considered as a vertex equivalent to the 
well-known edge decomposition (deletion-contraction relations). Edge decomposi- 
tion formulae of the form / (G) = a (e) / (G — e) + /? (e) / (G/e) apply to the Tutte 
polynomial and derived graph invariants, for instance the number of spanning trees 
or the rehability polynomial. Indeed, it was shown by J.G. Oxley and D.J. A. Welsh, 
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[41j , that the Tuttc polynomial is in a certain sense universal, meaning that all other 
graph invariants that satisfy edge decomposition formulae can be derived from the 
Tutte polynomial by substitution of variables. A textbook presentation is given in 
[lOj . A general framework analyzing universality properties of graph polynomials 
is studied in [25] . 

It seems natural to ask for the most general vertex decomposition formula. Let 
us assume that we try to construct an ordinary generating function / (G) that 
counts some type of vertex induced subgraphs with respect to the number of ver- 
tices. Which properties should such a function have? If the subgraphs in question 
are composed from subgraphs of the components then we can expect multiplica- 
tivity of / with respect to components of the graph. In order to assign the value 
/ (G) uniquely to a graph G by application of a decomposition formula as given in 
Theorem 1131 certain initial values for the null graph and the empty graph have to 
be given. Therefore, we presuppose the following four properties of /: 

(a) (Multiplicativity) If Gi and G2 are components of G then / (G) = / (Gi) / (G2). 

(b) (Recurrence relation) Let a, /3, 7 € M and let w be a vertex of G, then 

(2) / (G) = af{G-v)+Pf{G-N [v]) + 7/ {G/v) . 

(c) (Initial condition) There exists (5 G M such that f {%) = 5 for the null graph 
0= (0,0). 

(d) (Initial condition) There exists e G R such that f {Ei) = e for a graph 
El = {{v} , 0) consisting of one vertex. 

Furthermore, in order to make / a well-defined graph polynomial, the result of 
computing / has to be the same, irrespective of the order in which we apply the 
enabled computation steps. In particular, it has to be independent of the order of 
the vertices, which we use to apply the decomposition formula (b). In general we 
may choose a, (3, 7, 5, e from a field of characteristic zero or from a ring. A graph 
invariant is proper if there are two graphs Gi and G2 with the same number of 
vertices such that /(Gi) ^ /(G2). 

Applying the conditions (b), (c), and (d) we obtain from Ei — v = Ei — N [v] = 
Ei/v = the equation 

e ^ {a + (3 + j) 5. 

Computing / {E2) — f {K2} in two ways using (a) and (b), respectively, results in 

= {a + P + j)e. 
Consequently, the values of 6 and e are determined: 

S = 1 

e = a + P + j 

If the constants a, (3, 7 are properly defined then the value of / (G) does not depend 
on the choice of the vertex v in equation ((2]). Consequently, the function / (G) does 
not depend on the order of the vertex decomposition ([2|). The calculation of / (P3) 
for path of three vertices yields, in case we start from a vertex of degree 1, 

/(P3) = (a + 7)'(a + /3 + 7) + /3(a + 7) + /3(a + /? + 7)- 
If we begin the vertex decomposition at the vertex of degree 2 then we obtain 

/ (P3) = a (a + /3 + 7) V /3 + 7 (a + 7) (a + /3 + 7) + /57 
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These two results coincide if /? 0, a = 1, or a + /3 + 7 = 1. In any case, there 
remain only two variables that can be chosen independently. In case of /3 = 0, all 
graphs with the same number of vertices result in the same polynomial. Therefore, 
this case does not yield any interesting applications. Ifa + /3 + 7 = 1 then / (G) = 1 
for all graphs. The only remaining choice, a = 1, gives for P = x (y — 1) and j = x 
the subgraph component polynomial Q{G; x,y). 

From this we get that Q{G; x,y) is universal among polynomials recursively 
defined using vertex deletion, vertex extraction and vertex contraction. More pre- 
cisely, we have the following theorem. 

Theorem 21 (Universality of Q{G] x, y)). (1) For a graph polynomial f{G; a, (3, 7, 
to be proper and well-defined by conditions (a)-(d) we have a = 1, 5 = 1 
and e = 1 + /3 + 7. 

(2) There is a unique proper graph polynomial U{G; f3,j) which is well-defined 
by conditions (a)-(d) and we have 

Q{G-x,y) = U{G-x{y-\),x) 

and 

C/(G;/3,7) = Q(G;7,- + 1) 

7 

5.4. Vertex Eliminations vs Edge Elimination. The subgraph component 
polynomial Q{G;x,y) can be regarded as counting vertex set expansions. In the 
literature there is a variety of graph polynomials, including the Tutte polynomial, 
which can be defined by counting edge set expansions. 

We have seen in Theorem 12 II that Q(G;x,y) is universal among the polynomials 
defined recursively via deletion, extraction and contraction of vertices. In [H [5] the 
polynomial S,{G;x,y, z) was shown to be universal among the polynomials defined 
recursively via deletion, extraction and contraction of edges. In this section we 
will show the connection of G(G; x, y) to the universal edge elimination polynomial 
i{G]x,y,z). 

The polynomial S,{G;x,y, z) generalizes both the Tutte and the matching poly- 
nomials, as well as the bivariate chromatic polynomial of [18 . We shall use the 
recursive decomposition of ^{G;x,y, z) from [5]: 

^(G; X, y, z) ^ ^{G - e; x, y, z) -\- y^{G/e; x, y, z) + z£,{G f e; x, y, z) 
^(Gi U G2; x, y, z) = C(Gi; y, z)^(G2; x, y, z) 
^{Ei;x,y,z) = x 

(3) m) = 1 

where Gi U G2 denotes the disjoint union of graphs Gi and G2, and the three edge 
elimination operations are defined as follows: 

Edge deletion:: We denote by G — e the graph obtained from G by simply 

removing the edge e. 
Edge extraction:: We denote by G j e the graph induced by \ {u, u} pro- 
vided e = Note that this operation removes also all the edges 
adjacent to e. 

Edge contraction:: We denote by G/e the graph obtained from G by uni- 
fying the endpoints of e. 
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We will rewrite the decomposition of Q[G] x, y) using Theorem [T3l 

Q(G; X, y) = Q{G - v; x, y) + xQ{G/v- x, y) + x{y - l)Q{G - N[v]]x, y) 
Q{Gi U G2;x, y) = Q{Gi;x, y)Q{G2;x, y) 
Q{Ei;x,y) ^ xy + 1 

(4) Q(0) = 1 

Theorem 22. Let G = {V,E) be a graph. Let L{G) — (T4,£^e) denote the line 
graph of G. Then the following equation holds: 

^{G-l,x,x{y-l)) = Q[L{G)-x,y) 

Proof. First, let us analyze the correspondence of the edge elimination operations 
in a graph to the vertex elimination operations in its line graph. Let Ve G Ve be the 
vertex in the line graph that corresponds to the edge e G i? of the original graph. 
By the definition of the edge and vertex elimination operations: 

(5) L{G - e) = L{G) ~ v, 

(6) L{G/e) = L{G)/v, 

(7) L{GU)^L{G)-N[v,] 

First let us check the connected graphs with up to one edge: 

If Ge {0,Si}, L(G) = 0, 
The equivalence ^(G; l,x, x{y — 1)) = 1 = (5(0) holds. 

If G is a single point with a loop, or G = P2, L{G) is a singleton, The equivalence 
C(G; 1, X, x{y - 1)) = 1 + X + x{y - 1) ^ I + xy = Q{Ei) holds. 

Next, we note that L(GiUG2) = L(Gi)UL(G2). Therefore, if the theorem holds 
for graphs Gi and G2, then it holds also for Gi U G2. Finally, the theorem follows 
by induction on the number of edges, using the decomposition formulae (j4|) and ([3|) 
and the correspondence of edge and vertex elimination operations. □ 

Problem 23. How does the distinguishing power of $^(G; x,y, z) compare to the 
distinguishing power of Q{G; x, y) ? 

5.5. Subset Expansion and Definability in Logic. Q{G;x,y) was defined as a 
generating function. Let us rewrite the definition of Q{G; x, y) in a slightly different 
way. Instead of summation over the number of the used vertices i , and the number 
of induced connected components j, we shall summate over all the possible subsets 
of vertices: 

(8) Q(G;x,y)= 5]xl^ly'=(C[-41). 

A<ZV 

This is a subset expansion formula, a term coined in [42j . The relationship 
between recursive definitions of graph polynomials and the existence of subset ex- 
pansion formulas has been studied from a logical point of view in |25| . Subset 
expansion formulas can often be used to show that a graph polynomial is definable 
in Monadic Second Order Logic, as studied in [28l |3T] . However, the exponent 
k{G[A\) in Equation ^ causes a problem, to remedy this, we use, like in [29] . 
an auxiliary order -< over the vertices. We will denote by F{A) the subset of the 
smallest vertices in every respective connected component. 
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(9) Q{G;x,y)= E ( 11 ^) ( 11 2^1 

ACV \veA I \u(iF{A) ] 

Note that the resuh does not depend on the used auxiliary order. 

Without having to go in the details of graph polynomials definable in Monadic 
Second Order LogicQ, Equation ^ shows that Q(G\ x, y) is a graph polynomial 
definable in Monadic Second Order Logic for graphs G = (V, E) with universe V 
and a binary edge relation. Therefore all the theorems from [251 [31] can be applied. 
In particular, the Feferman-Vaught-type theorems from [28] guarantee existence of 
reduction formulas like multiplicativity from Theorem I12[ or the one in Theorem 
[18] for the join, not only for the disjoint union or the join operation, but for a wide 
class of MSOL-definable operations. Also, a general theorem from [32J guarantees 
the existence of recurrence formulas, as proven in Propositions [TB] and 1171 for a 
wide class of recursively defined families of graphs, as studied also in [3S|. Among 
these we have the wheels Wm the ladders L„ and the stars Stavn- It should not be 
difficult to compute the recurrence relations for these explicitly. 

We shall exploit MSOL-definability also for our complexity analysis in Section 
[92] 

6. Clique Separators 

The simplest case of a clique separator in a graph G is an articulation, i.e. a 
vertex whose removal from G results in an increase of the number of components. 
Let V be an articulation of G and let H and K be subgraphs of G such that 
G = HUK and H D K = {{v} , 0). It is well known, cf. [10], that in this case the 
Tutte polynomial T{G; x, y) satisfies 

T{G-x,y)^T{H-x,y)-T{H-x,y) 

In the case of the subgraph component polynomial the situation is a bit more 
complicated: 

Theorem 24. Let v be an articulation of G and let H and K he subgraphs of 
G such that G ^ H U K and H C\ K ~ ({u} ,0). Then the subgraph polynomial 
Q (G) = Q (G; X, y) satisfies 

Q{G) = Q{H -v)Q{K -v) 

+ — [Q (H) ~Q{H- v)] [Q (K) ~Q{K-v)]. 
xy 

Proof. The first product of the polynomial is the generating function for the number 
all vertex induced subgraphs that do not contain the articulation v. The product 
is justified by Theorem [T2l The second term counts all remaining subgraphs, i.e. 
those ones containing vertex v. Here the equation ([T]) from the proof of Theorem 
[T^ has to be modified. The vertex v is counted twice because it belongs to both K 
and H. This double counting is corrected by multiplication with x^^. By analogy, 
we introduce the factor y~^ in order to avoid that the component containing v is 
counted twice. □ 



^ The interested reader can consult |21| for the use Monadic Second Order Logic in finite model 
theory, and |13) for its use in graph theory. 
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Theorem [24] can be generalized in order to cover clique separators with more 
than one vertex. Let G = {V, E) be a connected graph and iJ, K subgraphs of 
G such that H n K = Kr and H U K = G. In this case Kr = {U, F) forms a 
separating clique of G. Here we assume that neither H nor K coincides with K^- 
The subgraphs H and K are called split components of G with respect to K^- 

Theorem 25. Let Kr = {U, F) be a clique separator of G such that there are 
two split components H and K . Then the subgraph polynomial Q (G) ~ Q (G; a;, y) 
satisfies 

Q{G) = Q{H -U)Q{K -U) 

+ ^ E i E E [~ir^''''Q{H-B)QiK-G). 

^ dl^AQU BDU\ACDU\A 

Proof. First we count all subgraphs that are induced by vertex subsets included 
in V \ U. These subgraphs are also subgraphs oi G — U. From Theorem [T^] we 
obtain Q [H — U)Q [K — U) as generating function for all subgraphs of G induced 
by subsets of V\U. 

For each subset A C U, let /y {H, A) be the number of vertex subsets X CV (H) 
of cardinality i with A C X such that the induced subgraph H [X] has exactly j 
components: 

{H,A)^\{X :ACX CV (H) A\X\ = tAk{H [X]) = j}\ 
The polynomial 

n n 

F{H,A)^Y.T.h^H^'^)'^'y' 

i=0 i=0 

is the ordinary generating function for the numbers fij {H, A). We define the num- 
bers fij {K, A) and the corresponding generating function F {K, A) for the second 
split component analogously. Let X C V (G) be a vertex subset with X nU = A. 
Then the component of G [X] that contains A is counted in F {H, A) and in 
F(K,A). There is indeed only one component counted twice, since A induces 
a clique of H and K, respectively. Thus we obtain 

(10) Q{G)=Q{H-U)Q{K-U) + - J2 -^F{H,A)F{K,A). 

The factor x"''^' takes into account that all vertices of A contribute to F{H,A) 
and to F {K, A). For each subset B C U, the subgraph polynomial oi H — B can 
be represented as a sum of generating functions as follows: 

Q{H-B)^ J2 F{H,A) 

ACU\B 

We define Q {H,U \B) ^ Q {H - B) and obtain 

QiH,U\B)^ J2 F{H,A) 

A<ZU\B 

or 

Q{H,B)^ F{H,A). 

A<ZB 
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By Mobius inversion, we obtain 

F{H,A) = ^ (-l)l-^H^Ig(ff,B) 

BCA 
BCA 
U\BCA 

^ ^ (-1)1^1 

Similarly, we can prove for each A C U that 

A) = (-1)1^1-1^1 ^ (-i)i^'g(i^-i3). 

The substitution of F {H, A) and F (K, A) in equation (HU]) yields 
Q(G) = Q{H -U)Q{K ~U) 



□ 



Problem 26. Can we have an analogue of Theorem \25\ for the case where the 
separating vertex set is not required to be a clique? 

7. Reconstruction 

The famous graph reconstruction conjecture by Kelly and Ulam 44J states the 
every undirected graph with at least three vertices can be reconstructed from a deck 
of its vertex-deleted subgraphs (more precisely from the corresponding isomorphism 
classes). See for example the papers [TH [TT] as an introduction into this field. 
Despite the fact that the conjecture is still open, many graph invariants and graph 
polynomials (e.g. the Tutte polynomial) are known to be reconstructible. We can 
show that also the subgraph component polynomial can be reconstructed from the 
deck of the subgraph polynomials of the vertex-deleted subgraphs. 

Theorem 27. The subgraph polynomial Q (G; x, y) for a graph G = (V, E) with 
n — \V (G)| > 3 is uniquely determined by the set of polynomials 

{Q{G~v-x,y):veV{G)}. 

Let Cj be the smallest power of y that appears at least twice among the terms x"~^y^ 
of the polynomials Q (G — v]x,y). Define 

n if uj = n — 1, 
uj else. 

Then the subgraph polynomial is given by 



Q{G;x,y)=x'' 



/%"-i^QfG-i;;i,y) dt 
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Proof. Let k (G) denote the number of components of a graph G. In each graph 
G — (y, E) with at least three vertices and at least one edge there exist two vertices 
u,v G V such that k{G — u) = k{G — v) = k (G). If the term x"~^y^ appears in 
the polynomial Q (G v; x, y) then the number of components of G — w equals j. 
Since k{G — v) > k (G) for each vertex v G V, the the smallest power of y that 
appears at least twice among the terms x^^^^y^ of the polynomials Q {G — v; x, y) is 
equal to k (G). There is only one exception: If G is the empty (edgeless) graph then 
the removal of each vertex of G decreases the number of components by one, which 
is taken into consideration within the definition of uj. Consequently, we obtain 
LJ = k{G). 

Each vertex induced subgraph with i < n vertices is counted exactly n ~ i times 
in the polynomial 

'{G-v;x,y). 



vev(G) 



The coefficient of f^y^ in 



G- 



equals i times the number of vertex induced subgraphs of G with exactly n — 
i vertices and j components. The integration with respect to t transforms t^~^ 
into ii* such that the vertex induced subgraphs are enumerated correctly by the 
coefficients of the resulting polynomial. Finally, the bounds of integration and 
the multiplication with performs the back-substitution in order to obtain an 
ordinary generating function with variables x and y. □ 

8. Random Subgraphs 

Now we assume that the vertices of G = {V,E) fail stochastic independently 
with a given (identical) probability q — 1 — p- We obtain the probability Pk (G) 
that a vertex induced subgraph of G has exactly k components from the subgraph 
polynomial: 



(11) 



y=o 



The sequence {Pk (G)}j,gj^ is the distribution of the number of components. Con- 
sequently, we obtain 

n 

Y.Pk{G) = l. 

fe=0 



-<3 



Figure 8. A 3 x 3 grid graph 



Figure [5] shows the distribution for the graph presented in Figure [51 
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Figure 9. Distribution of the number of components 

The probability Pi (G) is cahed the residual connectedness reliability. Boesch, 
Satyanarayana, and SufFel ^9] showed that the computation of Pi (G) is a #P- 
hard problem, even in planar bipartite graphs. Since Pi (G) can be obtained in 
polynomial time from the subgraph polynomial by applying the relation (jlip . we 
obtain the following statement. 

Corollary 28. The computation of the subgraph polynomial is a #P-hard problem. 
It remains #P-hard for the class of all planar bipartite graphs. 

9. Computational complexity of Q{G; x, y) 

9.1. Complexity of evaluation. We have already seen in Corollary[28lthat Q(G; x, y) 
is tJP-hard to compute. Now we deal with a problem of evaluation oi Q{—]x,y) at 
a given point (x, y) G for arbitrary input graph G. 

Theorem 29. For every point {x,y) G Q^, possibly except for the lines xy — 0, 
y ^ 1, X — —1 and x = —2, the evaluation of Q(G\x,y) for an input graph G is 
\\P-hard. 

C. Hoffmann in [26] showed the following: 

Theorem 30 (Hoffmann 2008). For every point {x,y,z) G , except possibly 
for the subsets x = 0, z — ~xy, {x,z) e {(1, 0), (2, 0)} and y £ {—2,-1,0}, the 
evaluation of ^{—;x,y, z) for an input graph G is \j,P-hard. 
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Proof of Theorem We use Theorem [30] and our Theorem ([22]) . Under the con- 
ditions of Theorem ([22]) . Hoffmann's exception sets are mapped to the hnes xy — Q, 
y = 1, X ^ —1 and x = —2. It follows that for every point {x,y) € that does 
not lay on one of those lines, the polynomial Q{—', x, y) is t|P-hard to evaluate even 
for an input line graph L{G). □ 

The evaluation of Q(—]x,y) is polynomial time computable for xy — and for 
y ~ I. It remains open whether it is polynomial time computable for a; = — 1 
and X — —2. One can also ask, whether there is some point (x^y) € Q^, in which 
Q(— ; X, y) is hard to evaluate for general input graph, but easy for input line graph. 

9.2. Parameterized complexity. Here we discuss the computational complex- 
ity of Q{G;x,y) for input graphs of bounded tree with, and for input graphs for 
bounded clique width. We do not need the exact definitions here. For background 
on tree- width the reader can consult [17]. Clique- width was defined in [T^. Both 
are discussed in [28] . 

Recall that the subgraph component polynomial is definable using the MSOL- 
formalism (definition [S]) with auxiliary order, while the result is order-independent. 
Hense, using a general theorem from [29l[28], we have 

Proposition 31. Q(G;x,y) is polynomial time computable on graphs of tree-width 
at most k where the exponent of the run time is independent of k. 

Moreover, applying the result of Courcelle, Makowsky and Rotics [T4], combined 
with the results from [40], we have a similar result for graphs of bouded clique 
width: 

Proposition 32. Q{G\x,y) is polynomial time computable on graphs of clique- 
width at most k where the exponent of the run time is independent of k. 

The drawback of the general methods of [29l [28] and [T^ , lies in the huge hidden 
constants, which make it practically unusable. However, an explicit dynamic algo- 
rithm for computing the polynomial Q{G; x, y) on graphs of bounded tree-width, 
given the tree decomposition of the graph, where the constants are simply expo- 
nential in k, can be constructed along the same ideas as presented in [43l[22]. For 
the graphs of bounded clique width, given the clique decomposition of the graph, 
we know an algorithm with constants doubly-exponential in A:. It is open whether 
an algorithm with constants simply exponential in k exists. For a comparison 
of the complexity of computing graph polynomials on graphs classes of bounded 
clique- width, cf. [33] . 

10. Conclusions and Open Problems 

We have shown that Q{G] x, y) is a universal vertex elimination polynomial. We 
have given a few combinatorial interpretations of its evaluations and coefhcients. 
We have proven various splitting formulas for Q(G;x,y) such as the multiplicativ- 
ity. Theorem [151 and Theorem]^!] Problem [^ asks for more such theorems. Besides 
having algorithmic importance, such splitting formulas increase our structural un- 
derstanding of the graph polynomial under study, and may help us in analizing its 
distictive power. 

We have looked at the graph polynomial Q{G; x, y) from various angles and com- 
pared its behaviour and distinguishing power with the characteristic polynomial. 
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the matching polynomial the Tutte polynomial and the universal edge elimination 
polynomial. We have not discussed the relationship of Q{G;x,y) to other graph 
polynomials, such as the interlace polynomial, [21 H], or the many other graph 
polynomials listed in [31j . 

We have seen that Q{G; x, y) distinguishes between graphs where these polyno- 
mials do not. We have not found cases where these other polynomials do distinguish 
between graphs where Q(G; x, y) does not. This is probably due to our lack of com- 
puterized tools for searching for such cases, cf. Problem [51 In Problem [221 we ask 
about comparing distinguishing power of Q{G]x,y) and the universal edge elimi- 
nation polynomial ^(G; x, y, z). This seems to be more tricky. We have given a few 
examples of graphs and graph families which are determined by Q{G; x,y). 

Problem 33. Find more graph invariants which are determined by Q(G;x,y) . 

Problem 34. Find more classes of graphs which are determined by Q{G;x,y). 

Returning to our motivation, we have only studied the simplest case of commu- 
nity structure in networks. We have studied the generating function of induced 
subgraphs with i vertices which have j components. More generally, one would 
want to study community structures where components are replaced by maximal 
fc-connected components. 

Problem 35. What are the appropriate generating functions which capture the 
essence of various community structures? 
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